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Ab initio calculations for the giant magnetoresistance (GMR) in Co/Cu, Fe/Cr, and Fe/Au mul- 
tilayers are presented. The electronic structure of the multilayers and the scattering potentials of 
point defects therein are calculated self-consistently. Residual resistivities are obtained by solving 
the quasi-classical Boltzmann equation including the electronic structure of the layered system, the 
anisotropic scattering cross sections derived by a Green's function method and the vertex correc- 
tions. Furthermore, the influence of scattering centers at the interfaces and within the metallic layers 
is incorporated by averaging the scattering cross sections of different impurities at various sites. An 
excellent agreement of experimental and theoretical results concerning the general trend of GMR 
in Co/Cu systems depending on the type and the position of impurities is obtained. Due to the 
quantum confinement in magnetic multilayers GMR can be tailored as a function of the impurity 
position. In Co/Cu and Fe/Au systems impurities in the magnetic layer lead to high GMR values, 
whereas in Fe/Cr systems defects at the interfaces are most efficient to increase GMR. 

PACS numbers: 73.63.-b,75.70.Cn,75.30.Et,71.15.Mb 
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A large number of experimental and theoretical inves- 
tigations was initiated by the discovery of giant magne- 
toresistance (GMR) in magnetic multilayers 0,0 to elu- 
cidate the microscopic origin of the phenomenon. Several 
authors have shown |3J, |j] that GMR in magnetic multi- 
layers is strongly influenced by changes in the electronic 
structure, especially the Fermi velocities, of the system 
in dependence on the relative magnetization alignment 
in adjacent layers. In realistic samples, however, spin- 
dependent scattering is considered to cause GMR. The 
effort to tailor GMR systems with high rates was accom- 
panied by a variety of experiments [3, la] and calculations 
[21 la, Loll which investigated the influence of dusting and 
doping by impurities. Agreement was reached concerning 
the dominant role of interface scattering || . The results 
of Marrows et al. Q, however, demonstrated the strong 
dependence of GMR on the position of the impurities 
with respect to the interfaces and on the valence differ- 
ence between impurity and host. The aim of this paper 
is to present ab initio calculations for the scattering cross 
sections and resulting GMR ratios in dependence on de- 
fect material and position for the standard systems of 
magnetoelectronics Co/Cu and Fe/Cr. In addition, the 
system Fe/Au was investigated because the system shows 
a high interface quality 10] and the structural properties 
should be closest in experiment and theory. The influ- 
ence of defects on interlayer exchange coupling in Fe/Au 
was investigated earlier jl l| . 



All calculations are performed within the framework 
of density functional theory in local spin density approx- 
imation applying the Screened KKR (Korringa-Kohn- 
Rostoker) Green's function method 0|. We have chosen 
a multilayer geometry in the so-called first antiferromag- 
netic (AF) maximum of interlayer exchange coupling for 
the Co/Cu system consisting of 9 monolayers (ML) Co 
separated by 7 ML Cu, denoted as C09CU7. The struc- 
ture of the superlattice was assumed to be an fee lattice 
with a lattice constant of 6.76 a.u. grown in (001) direc- 
tion. Each atomic plane is represented by one atom in 
the prolonged unit cell with 32 atoms. A similar config- 
uration was chosen for Fe / Cr and a perfect bec stacking 
in (001) direction with a lattice constant of 5. 50a. it. was 
assumed. This is samewhat larger than the lattice con- 
stant of Fe and Cr, but provides the correct magnetic 
order within the LDA and the atomic sphere approxima- 
tion for the potentials The Fe/Cr system consists of 
9 ML Fe and 9 ML Cr which possesses an AF g round state 
confirmed by calculation and experiment [lfj . Structural 
relaxations at the interfaces were neglected. Despite of 
the AF order the Cr layer will be refered to as the non- 
magnetic layer in the system, to use the same notation 
as for the other system under consideration. We consider 
Fe/Au multilayers with the same structural data reported 
by the experiments 0] that is a structural bcc-fcc transi- 
tion. The lattice constant for Fe is ab cc — 5. 4163a. it., and 
for Au a f cc — \/2ai, cc . The thickness of 9 ML of the Au 
layer was chosen in accordance with the experimentally 
obtained interlayer exchange coupling strength which fa- 
vors an antiparallel coupling for this Au thickness. 

The self-consistent electronic structure of the ideal 
host, without any impurities, is described by the one- 
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particle Green's-function, whose structural part G^t'lM) 
is expanded into a site and angular- momentum basis 

The new aspect for superlattices is that we consider 
a lattice with a basis. The index n is now a shorthand 
notation for lattice vector Rat and basis vector of the 
atoms in the unit cell. To simulate a substitutional point 
defect one atom in the lattice is replaced by another. The 
site index of the impurity position is denoted by /i. 

The impurity Green's function G™y (E) of the multi- 
layer including an impurity atom at a defined position is 
obtained by the solution of an algebraic Dyson equation 
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The ^''-summation involves all sites in the vicinity of the 
site fi for which the differences of the single site t-matrices 
Atjy/ = t\,, — 1\„ of the multilayer with and without de- 
fect are significant. The single site t-matrices are derived 
from the angular-momentum dependent scattering phase 
shifts of the potentials in atomic sphere approximation 
(ASA). The differences Ai^ characterize the potential 
perturbation caused by the defect. In the calculations 
we take into account angular momenta lmax < 3. We al- 
low for potential perturbations up to the second atomic 
shell around the impurity atom. Charge multipoles up to 
lmax — 6 are taken into account. Since the systems un- 
der consideration are magnetic all properties mentioned 
above depend also on spin quantum numbers a =] , J, for 
majority and minority electrons, respectively. 

Using the impurity Green function, the self- 
consistently calculated potential perturbation and the 
wave function coefficients of the superlattice Bloch states 
we derive the microscopic spin-conserving transition 
probability P kk , for elastic, that is on-shell-scattering of 
a Bloch wave k into a perturbed Bloch wave k'. kis now 
a shorthand notation for the wave vector k and band in- 
dex v, a denotes the spin quantum number [lfiLlltij . The 
transition probability is given by Fermi's golden rule 
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The formalism is restricted to dilute alloys since we as- 
sume a linear dependence with the number of defects 
cN. Furthermore, spin-orbit coupling and the resulting 
spin-flip processes are neglected by the non-relativistic 
scheme. 

The transition matrix elements T kk , for the scattering 
of Bloch electrons by an impurity cluster embedded in 
an otherwise ideal translational invariant multilayer are 
given by [3, E3| 



T kk> = \ E Cl{k,a)T^Ci{k',cj). 



(3) 
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C£(fc,cr) are the expansion coefficients for the superlat- 
tice wave function in an angular momentum basis. V 



denotes the total crystal volume. Using spherical poten- 
tials (ASA) the matrix elements T£2> are derived from 
the structural Green's function matrix elements G™ L , of 
the perturbed system and the potential perturbation At? 

mum- 

Summation over all final states leads to the spin and 
state dependent relaxation time 
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which depends on spin u, Bloch state k and im- 
purity position r M in the superlattice. The k and 
spin-dependence of the scattering rates is treated fully 
quantum-mechanically without adjustable parameter. 
The dependence on the effect position is involved by 
the implicit dependence of the impurity Green's function 
G^l' on the r M . Up to this point we consider a dilute 
alloy of impurity atoms all of them occupying a chosen 
site fi in the unit cell. That is, the alloying is restricted 
to certain atomic planes in the multilayer. These planes 
correspond to the impurity <5-layers experimentally inves- 
tigated by Marrows and Hickey Q. 

The conductivity is calculated by solving the quasi- 
classical Boltzmann equation EH- Thus, the vector mean 
free paths are obtained by 
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This includes besides to the anisotropic relaxation times 
as the second term on the r.h.s the computational de- 
manding scattering-in term (vertex corrections) which 
completes the description of impurity scattering |lfjl2(l]. 
The band structure is included via Fermi velocities v k 
and the k' summation over all states on the anisotropic 
Fermi surface. The impurity scattering enters via the 
relaxation times r k (fi) and the microscopic transition 
probabilities P kk , (/x) . Eq. (0 is solved iteratively to de- 
termine the vector mean free path A£(/i) of an electron 
with spin a in a state k. To our knowledge, up to now 
the semi-classical calculations of the GMR have been 
mostly performed within relaxation time approximation 
only 0,0,0,113, 001 , which neglects the scattering-in 
term. Zhang and Butler proposed a simplified method to 
include the vertex corrections by renormalization of the 
electron life times using an adjustable parameter [2fi| . 
The deviation of the relaxation time approximation in 
comparison with the solution including the vertex cor- 
rections are discussed in Fig. [2 

Based on the solution of Eq. © the spin-dependent 
conductivity tensor a° ' [p) is given by a Fermi surface 
integral [lg] 
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PIG. 1: Conductivity of Cog/Cu?, Fe<jjCri ) and Fecj/Auj for P and AP alignment in dependence on the position of 'self 
impurities, respectively; a,c,e) Assuming scattering at the inserted 5-layer only, b,d,f) assuming 5-layer scattering (50%) and 
interface scattering (50%) 



With Mott's two current model [26] the total conductiv- 
ity cr(/i) = cr T (/i) + &*(fj,) is obtained by spin summa- 
tion. The CIP resistivity is obtained by the inverse of 
the CIP conductivity due to the diagonal structure of 
the conductivity tensor for the tetragonal systems under 
consideration 



fft(^) + trl(/i) 



(7) 



To describe the existence of an overall distribution of 
impurities in the multilayer the transition probabilities 
of the different <5-layers have to be superimposed. Fol- 
lowing this idea, layer-dependent relaxation times t£ (p) 
are added 
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including weighting factors x((i) that account for the rel- 
ative concentration of defects at the corresponding posi- 
tions [i in the unit cell. 

The most driving aspect of magnetic multilayers is the 
drastic change of the conductivity cr as a function of the 
relative orientation of the magnetic layer moments, par- 
allel (P) or anti-parallel (AP). The relative change defines 
the GMR ratio 



GMR[p) = 
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CONDUCTIVITY AND QUANTUM WELL 
STATES 



We consider the multilayers described above and in- 
vestigate the scattering properties of impurities of the 
nonmagnetic component in the magnetic layers and vice 
versa. 

The analysis of the transport coefficients is focused on 
the current-in-plane (CIP) geometry. The total CIP- 
conductivities normalized to the defect concentration c 
caused by impurities of the magnetic material at differ- 
ent positions in the nonmagnetic layer and nonmagnetic 




Sc Ti V CrMnFeCoNiCuZn 

FIG. 2: Influence of vertex corrections: Dependence of the 
spin dependent resistivity for Cog/Cur multilayers on 3d im- 
purities at the interface position of the magnetic layer p|| cal- 
culated including vertex corrections and p\\ without 



defects in the magnetic layer are shown for all considered 
systems in Fig.^,,c,e), respectively, for both configura- 
tions of the magnetization directions (P, AP). The con- 
ductivity differs by orders of magnitude as a function of 
impurity position (keep in mind the logarithmic scale). 
The largest values occur for impurity positions where the 
quantum confinement produces many Bloch states with 
low probability amplitude. The eigenstates show strong 
quantum confinement due to the superlattice potential. 
That is, the probability amplitude is modulated by the 
layered structure and can even tend to zero at partic- 
ular sites of the supercell 0, 0|- The consideration of 
additional interface scattering in Fig.Q]b,d,f) will be dis- 
cussed below. 

The influence of vertex corrections is quantified in 
Fig. [21 For different defects at the Co interface position 
in Co/Cu the deviation of the resistivity p\\ calculated 
including vertex corrections in Eq. (2J and without (pu) 
are shown for both spin bands. In the minority channel 
(upwards triangles) the deviations are of the order of few 
percent (except for Cr). In the majority channel the ver- 
tex corrections are large for defects where sp-scattering 
dominates. This is the case for impurities with a similar 
electronic structure like the host, at least for this spin 
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direction. For Sc, Ti, V, and Cr with an opposite mag- 
netic moment with respect to Co and a dominating sd 
scattering the vertex corrections are small as for the mi- 
nority spin direction. Summarizing, one can state, that 
the variations of the resistivity by incorporation of the 
vertex corrections in the solution of the Boltzmann equa- 
tion are less than 20% and do not change the qualitative 
behaviour of GMR. The neglect of vertex corrections may 
change the results quantitatively, but the general trend 
is conserved. 

To analyze the influence of the quantum confinement 
on the character of the eigenstates a classification of the 
eigenstates according to Ref. Q was performed. The dis- 
tribution of the probability amplitude in the different re- 
gions of the multilayer was analyzed and the states are 
labeld by the region with the highest averaged probabil- 
ity amplitude. The eigenstates are classified in 4 types 
and are labeled by quantum well states (QWS) in the 
magnetic layer Cm (Co, Fe), QWS in the nonmagnetic 
layer Cjv (Cu, Cr, Au), interface states Cj, and extended 
states Ce which have a probability amplitude of approx- 
imately the same size in all regions. To emphasize the 
contribution of the different classes of eigenstates the to- 
tal conductivity (Eq. was projected to the layers v 
in the supercell according to Ref. 
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with c1{v) the expansion coefficient of the Bloch eigen- 
state (k, a) at the site v in the supercell. The nor- 
malization of the Bloch state to the unit cell yields 
'Yl,v \ c ki l ')\ 2 ~ ^- I* 1 addition, the site dependent con- 
ductivity a a (y) was split into the contributions of the 4 
typical classes of eigentstates which are shown in Fig. 
It is evident, that a large contribution of the CIP con- 
ductivity is carried by the QWS in the magnetic and 
nonmagnetic layer. Especially, in the minority channel of 
the Co/Cu system and the majority channel of the Fe/Cr 
system the conductivity is dominated by contributions of 
the magnetic quantum well states. By considering only 
a few types of defects in the sample the conductivity is 
dominated by one spin direction in most cases. Assuming 
more types of defects and taking into account the effect 
of self averaging this tendency is reduced. 

Due to the quantum size effects the relaxation times 
show a strong variation for the states at the Fermi level 
which determine the conductivity. All Bloch states with 
a nearly zero probability amplitude at the impurity site 
undergo a weak scattering and cause extremely large re- 
laxation times. The state dependent relaxation times are 
distributed over several orders of magnitude, especially 
for defects inside the metallic layers (see Fig. 31 . This is 
a new effect that does not occur in bulk systems. The 
panels show the relative amount of relaxation times Tfc for 
the states at the Fermi level for the Co/Cu system (top 
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FIG. 3: Spin-dependent, layer-projected conductivity of 
CogCu7, FegCrr, and FegAu-? for P alignment with inter- 
face defects, typical classes are marked by grey - extended, 
black - magnetic QWS (Co, Fe), white - nonmagnetic QWS 
(Cu, Cr, Au), and light grey interface states. 



viewgraph) and the Fe/Cr system (bottom panel). The 
spin resolved histograms for Cu defects in bulk Co are 
given for comparison (topmost subpanel). The remain- 
ing subpanels represent relaxation times for Cu defects 
in the center of the Co layer and for Cu defects at the 
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Co/Cu interface. The color of the bars labels the char- 
acter of the eigenstates: extended multilayer states are 
shown in dark grey, QWS confined to the magnetic layer 
are shown in black. QWS confined to the nonmagnetic 
layer are given in white, and interface states are given 
by light grey bars. For defects inside the magnetic layer 
the maximum of the distribution coincides with that in 
the bulk material. In addition, a long tail for high values 



occurs caused by states which have a small probability 
amplitude at the defect position, e.g. quantum well states 
in the nonmagnetic layer or interface states. This is best 
seen in the middle subpanel for Cr defects in the center 
of the Fe layer for the Fe/Cr multilayer. Quantum well 
states confined to the Fe layer have a large probability 
amplitude at the defect position and as a result smaller 
relaxation times than in the bulk system. Interface states 
with a tail penetrating the Fe layer are scattered on an 
intermediate level and Cr quantum well states with the 
lowest probability amplitude at the defect position are 
scattered weakly. 

The states with large relaxation times although not nu- 
merous are highly conducting and nearly provoke a short 
circuit. This is the case for Co impurities in the Cu layer 
for the Co/Cu multilayer, compare Fig.Q] top panel, and 
for Fe defects in the Au layer for the Fe/Au system. This 
effect is mainly obtained for impurities in the center of 
the layers and is related to the fact that in-plane trans- 
port is mostly driven by quantum well states yj. This 
peculiar behavior of conductivity is in agreement with 
the results of Blaas et al. who found higher resistivi- 
ties for Co/Cu multilayers with interdiffusion restricted 
to the interface layers than for alloying with Cu atoms in 
the Co layers. 

For comparison with experiments we have to mention 
that the large absolute values would hardly be obtained 
experimentally since they correspond to idealized sam- 
ples with perfectly flat interfaces and defects at well de- 
fined positions in the superlattice. As soon as we con- 
sider an overall distribution of defects in the multilayer 
the highly conducting channels are suppressed. The gen- 
eral trend, however, survives. This phenomenon of highly 
conducting electrons confined to one layer of a multilayer 
structure is called electron wave guide or channe ling ef- 
fect |24l and was also experimentally verified |2flj . 

Structural investigations of Co/Cu multilayers on an 
atomic scale gave evidence that most of the struc- 

tural imperfections appear next to the interfaces. To in- 
vestigate the influence of more than one type of scattering 
centers in one sample a simplified defect distribution was 
assumed. In addition to the specific <5-layer defects of 
the magnetic layer material in the nonmagnetic interface 
atomic layer and defects from the nonmagnetic material 
in the magnetic interface layer are considered to simulate 
an intermixed region at the interface. For the concentra- 
tion weights x{p) entering Eq. © we choose 25% for de- 
fects in both of the interface layer and 50% for defects in 
the 5-layer and the resulting conductivities are shown in 
Fig.^b),(d),(f), respectively. First, an overall reduction 
of the resistivities is obtained, caused by the distribution 
of defects at different positions and the resulting higher 
probability that also quantum well states are scattered at 
one or the other type of defect. According to Eq. © the 
defect with the highest scattering rate l/r^ dominates 
the total relaxation time Tfc. 
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FIG. 5: Dependence of GMR in Cog/Cu?, Feg/Crr, and Feg/ 'Au^on the position of 'self impurities, a-c) Assuming scattering 
at the inserted 5-layer only, d-f) assuming 5-layer scattering (50%) and interface scattering (50%). The thin dashed lines marks 
the GMR ratio assuming interface impurities only, and the thick dashed line gives the GMR assuming isotropic relaxation 
times, the inset in panel (d) shows the experimental result from Ref. [f| 



GIANT MAGNETO-RESISTANCE 

The GMR ratios derived from the conductivities in 
Fig- Hare shown in Fig. 03 Assuming scattering centers 
in the <5-layer only (Fig. (a),(c),(e)) huge GMR ratios 
are obtained especially for Cu defects in the Co layer. 
Introducing additional interface scattering with a weight 
of 50% causes strongly reduced values (Fig. |^b),(d),(f). 
The thin dashed line in the lower panels of Fig. is the 
GMR ratio caused by interface scattering. This value 
would correspond to the reference value in the experi- 
ments of Marrows and Hickey without <5-layer [6]. The 
thick dashed line in Fig. Ogives the GMR value obtained 
with the assumption of a constant relaxation time with- 
out any spin or state dependence (intrinsic GMR). In 
comparison to this case of isotropic scattering the inser- 
tion of an additional (5-layer increases GMR, mostly at 
the interfaces. 

Comparing the trend of GMR an excellent agreement 
with the experiment (inset in Fig. is obtained for the 
Co/Cu system. To our knowledge, up to now similar ex- 
perimental investigations of the GMR dependence on the 
defects position are not available for Fe/Cr multilayers. 
The importance of the interface scattering to obtain a 
large GMR effect in Fe/Cr was pointed out by several 
authors |H & The weighting factors for the dif- 
ferent scattering mechanisms are derived from the struc- 
tural investigations by Davies et al. We choose 0.45 
for Cr defects in the Fe interface atomic layers and 0.05 
for Fe defects in the Cr interface layer. The contribution 



from the (5-layer is fixed to 0.5 as in the case of the Co/Cu 
multilayer. The most striking feature in comparison to 
Co/Cu is the reduction of the GMR ratio by introducing 
Cr defects in the Fe layers as shown in Fig. EJd) and e) . 

We still have to mention, that the calculated values 
are two orders of magnitude larger than the experimen- 
tal ones. The reason is the restriction to substitutional 
point defects. In addition to these much more scat- 
tering mechanisms are active in real samples. Assum- 
ing self-averaging the results could be corrected towards 
the experimental ones by an additional spin- and state- 
independent relaxation time r (thick dashed line in fig. 
2d-f) 0. 

In contrast to Ref. j^Hl the present results were ob- 
tained assuming the above described impurity distribu- 
tion only and are focussed on the impurity scattering 
rates only^Another difference to the experimental setup 
in Ref. [3j| is the considered geometry. The experi- 
mentally investigated samples have been Co/Cu/ Co spin 
valves grown on a buffer layer and protected by a cap 
layer. As a consequence the GMR ratios are nearly sym- 
metric as a function of the impurity position in the Cu 
layer but asymmetric for defects in the Co layer. The 
calculations are performed in supercell geometry which 
is reflected in the symmetry of the results with respect to 
the defect position in both layers, Cu and Co. A possi- 
ble influence of superlattice effects in metallic multilayers 
was shown to be negligible [3t| . 

Fig. El compiles the trend of GMR caused by 3d- 
,4s- (left column) and 4d-,5s-transition metal impurities 
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perturbations characterized by a spin-dependent scatter- 
ing strength t° ' . Comparing with these results one can 
state the following. The reader should note that the or- 
der of Co and Cu in Fig. 2 in ref. 0| is reversed in com- 
parison to Fig. UJi) and d) of this work. The insertion of a 
(5-layer in the Co layer enhances the GMR in comparison 
to the undoped case in both models. The differences of 
the results are caused by the approximation used for the 
scattering strength. In this work the scattering poten- 
tials were determined self consistently also for a region 
around the defects. Evaluating Eq. (jSJ the Born series 
expansion of the scattering operator T in terms of the 
single site t-matrix and all multiple scattering contribu- 
tions are included completely. 

Comparing the influence of 3d defects at the interface 
with that caused by an ordered interface alloy [38] a sim- 
ilar trend for the GMR is obtained. One should compare 
Fig. 5 (left column) in ref. j^| and the GMR values for 
the position x — of the (5-layer in Fig. El For the lighter 
3d elements up to Mn the GMR is lowered with the in- 
terface alloy, whereas for heavier elements the GMR is 
maintained or even increased. 

All the calculations are carried out for the limit of low 
defect concentration. That means the scattering at dif- 
ferent defects is treated independently. For a typical con- 
centration of 1% we obtain an imaginary part of the self 
energy of the order of 10 _3 eF . This is small in compar- 
ison with the typical band width in the transition met- 
als. In contrast to the work of Tsymbal et al. the 
contributions to the conductivity arising from interband 
transitions are expected to be small. 



SUMMARY 



FIG. 6: Dependence of the GMR on the impurity 5-layer 
doped with different transition and sp-metals and the position 
in the in Cog/Cu7(to-p panel), Feg/Cr7 (middle panel), and 
F eg I Am (lower panel); the label of each subpanel denotes the 
layer where the (5-layer is inserted; closed symbols: positions 
in the center, and open symbols mark the interface positions 



(right column) as a function of position in the magnetic 
layer. The CIP-GMR ratios are given for the correspond- 
ing defects in the middle of the magnetic and nonmag- 
netic layer, respectively, marked by the closed symbols. 
The open symbols correspond to a position of the (5-layer 
at the interface. The horizontal dashed lines indicate the 
case of interface scattering only. This means we consider 
a (5-layer of the magnetic material in the nonmagnetic 
interface atomic layer and vice versa. The thick dashed 
lines give the results assuming a spin and state indepen- 
dent relaxation time. 

In a previous work ,1| we used a simplified model to de- 
scribe the scattering. We assumed (5-peak like potential 



In conclusion, the self-consistent calculation of the 
scattering properties and the improved treatment of the 
Boltzmann transport equation including vertex correc- 
tions provide a powerful tool for a comprehensive theo- 
retical description and a helpful insight into the micro- 
scopic processes of CIP-GMR. The experimentally found 
trends concerning the doping with various materials at 
different positions in the magnetic multilayer could be 
well reproduced which means that spin-dependent impu- 
rity scattering is the most important source of GMR. The 
theoretical results show furthermore that interface scat- 
tering caused by intermixing plays a crucial role and has 
to be taken into account in any system under considera- 
tion. Selective doping of the multilayer with impurities in 
specific positions causes variations of GMR which could 
be well understood by the modulation of spin-dependent 
scattering due to quantum confinement in the layered 
system and by the spin anisotropy a. 
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